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I summarize the present status of global analyses of neutrino oscillations, including the most recent KamLAND 
and K2K data, as well as the latest solar and atmospheric neutrino fluxes. I give the allowed ranges of the three- 
flavour oscillation parameters from the current worlds' global neutrino data sample, their best fit values and 
discuss the small parameters a = Am SOL /Am ATM and sin #13, which characterize the strength of CP violation 
in neutrino oscillations. I briefly discuss neutrinoless double beta decay and the LSND neutrino oscillation hint, 
as well as the robustness of the neutrino oscillation results in the presence of non-standard physics. 



1. INTRODUCTION 

The discovery of neutrino oscillations by com- 
bining data from solar P, atmospheric [2], reac- 
tor [5j and accelerator 0] neutrino studies has 
marked a turning point in our understanding 
of nature and has brought neutrino physics to 
the center of attention of the particle, nuclear 
and astrophysics communities. This culminates 
a heroic effort which started over four decades 
ago and now firmly establishes the existence of 
small neutrino masses, confirming theoretical ex- 
pectations which date back to the early eight- 
ies. Neutrino mass arise from the dimension-five 
operator U4>4> where tf> the SU(2) <g> U(l) Higgs 
doublet and I is a lepton doublet jSJ. Nothing 
is known from first principles about the mech- 
anism that induces this operator, its associated 
mass scale or flavour structure. Its most pop- 
ular realization is the seesaw mechanism which 
induces small neutrino masses from the exchange 
of heavy states. Although inspired by unifica- 
tion, the effective model-independent description 
of the seesaw at low-energies is in terms of the 
SU{2) (g) U(l) gauge structure and contains a 
small effective Higgs triplet contribution to the 
neutrino mass (i . Such general seesaw has 24 
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parameters, 3 more than the tripletless case [H] 2 . 

The structure of the three-flavour lepton mix- 
ing matrix in various gauge theories of neutrino 
mass was given in .64 and it was also argued that, 
on general grounds, massive neutrinos should be 
Majorana particles, leading to neutrinoless dou- 
ble beta decay (P(3qu ) @- Current neutrino os- 
cillation data are well described by the simplest 
unitary lepton mixing matrix neglecting CP vio- 
lation. The effect of Dirac CP phases in oscilla- 
tions and Majorana phases in /?/3oi/ constitute the 
main challenge for the years to come. 

Here I focus mainly on the determination of 
neutrino mass and mixing parameters in neu- 
trino oscillation studies, a currently thriving in- 
dustry with many new experiments currently 
underway or planned. The interpretation of the 
data requires good calculations of solar and at- 
mospheric neutrino fluxes neutrino cross 
sections and experimental response functions, as 
well as a careful description of neutrino propa- 
gation both in the Sun and the Earth, includ- 
ing matter effects |12I13| . After summarizing the 
latest (post-Neutrino-2004) global analysis of 3- 

2 In the mass basis these correspond to the 12 mixing an- 
gles and 12 CP phases (both Dirac and Majorana-type) 
that characterize the full 3x6 charged current seesaw lep- 
ton mixing matrix ; 6j. These are exactly the same param- 
eters involved in the description of leptogenesis 7 , though 
in this case the use of the weak basis seems more conve- 
nient. Note also that current nomenclature (type-I versus 
II seesaw) is opposite from the original one in [H]. 
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neutrino oscillation parameters |T1| JS] I briefly 
discuss their impact on future /3f3o„ searches ^Sj • 
This writeup combines my plenary talk and the 
WG1 parallel session talk on non-standard sce- 
narios of neutrino conversion. 

2. TWO-NEUTRINO OSCILLATIONS 

Let us first consider neutrino oscillations in the 
two-flavour approximation |14j . 

2.1. Solar + KamLAND 

The solar neutrino data includes the measured 
rates of the chlorine experiment at the Homes- 
take mine (2.56 ± 0.16 ± 0.16 SNU), the most 
up-to-date results of the gallium experiments 
SAGE (66.9 ±11 ±|1 SNU) and GALLEX/GNO 
(69.3 ± 4.1 ± 3.6 SNU), as well as the 1496-day 
Super-K data in the form of 44 bins (8 energy 
bins, 6 of which are further divided into 7 zenith 
angle bins). The SNO data include the most re- 
cent data from the salt phase in the form of the 
neutral current (NC), charged current (CC) and 
elastic scattering (ES) fluxes, as well as the 2002 
spectral day/night data (17 energy bins for each 
day and night period). 

The analysis methods are described in \Fij\ an d 
references therein [T2|- We use a generalization 
of the pull approach for the \ 2 calculation origi- 
nally suggested in Ref. [20] m which all system- 
atic uncertainties such as those of the eight solar 
neutrino fluxes are included by introducing new 
parameters in the fit and adding a penalty func- 
tion to the x 2 ■ The method |14ll9j is exact to all 
orders in the pulls and covers the case of corre- 
lated statistical errors [21] as necessary to treat 
the SNO-salt experiment. This is particularly in- 
teresting as it allows us to include the Standard 
Solar Model 8 B flux prediction as well as the SNO 
NC measurement on the same footing, without 
pre-selecting a particular value, as implied by ex- 
panding around the predicted value: the fit itself 
chooses the best compromise between the SNO 
NC data and the SSM prediction. 

KamLAND detects reactor anti-neutrinos at 
the Kamiokande site by the process v e + p — > 
e + + n, where the delayed coincidence of the 
prompt energy from the positron and a charac- 



teristic gamma from the neutron capture allows 
an efficient reduction of backgrounds. Most of 
the incident v e flux comes from nuclear plants at 
distances of 80 — 350 km from the detector, far 
enough to probe the LMA solution of the solar 
neutrino problem. The neutrino energy is related 
to the prompt energy by E v — E pr + A — m e , 
where A is the neutron-proton mass difference 
and m e is the positron mass. For lower energies 
there is a relevant contribution from geo-neutrino 
events to the signal [221 - To avoid large uncer- 
tainties associated with the geo-neutrino flux an 
energy cut at 2.6 MeV prompt energy is applied 
for the oscillation analysis. 

First KamLAND results were taken from 
March to October 2002. Data corresponding to 
a 162 ton-year exposure gave 54 anti-neutrino 
events in the final sample, after all cuts, while 
86.8 ± 5.6 events are predicted for no oscillations 
with 0.95 ± 0.99 background events. The prob- 
ability that the KamLAND result is consistent 
with the no-disappearance hypothesis is less than 
0.05%. This gave the first evidence for the disap- 
pearance of neutrinos traveling to a detector from 
a power reactor and the first terrestrial confirma- 
tion of the solar neutrino anomaly. 

New KamLAND data were presented at Neu- 
trino 2004 [JjJ. With a somewhat larger fiducial 
volume of the detector an exposure correspond- 
ing to 766.3 ton-year has been obtained between 
March 2002 and January 2004 (including a re- 
analysis of the previous 2002 data). In total 258 
events have been observed, versus 356.2 ±23.7 re- 
actor neutrino events expected in the case of no 
disappearance and 7.5 ±1.3 background events. 
This leads to a confidence level of 99.995% for v e 
disappearance, and the averaged survival proba- 
bility is 0.686 ± 0.044(stat) ± 0.045(syst). More- 
over evidence for spectral distortion consistent 
with oscillations is obtained £|. 

It is convenient to treat the latest KamLAND 
data binning them equally in 1/E pr , instead of 
the traditional bins of equal size in E pv . Vari- 
ous systematic errors associated to the neutrino 
fluxes, backgrounds, reactor fuel composition and 
individual reactor powers, small matter effects, 
and improved v e flux parameterization are in- 
cluded [T2|. O ne finds a beautiful agreement be- 
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twecn KamLAND data and the region implied by 
the LMA solution to the solar neutrino problem, 
which in this way has been singled out as the only 
viable one in contrast to the previous "zoo" of os- 
cillation solutions |19I23| . However the stronger 
evidence for spectral distortion in the recent data 
leads to improved information on AmgQ L , sub- 
stantially reducing the allowed region of oscilla- 
tion parameters. From this point of view Kam- 
LAND has played a fundamental role in the res- 
olution of the solar neutrino problem. 

Assuming CPT invariance one can directly 
compare the information obtained from solar neu- 
trino experiments with the KamLAND reactor re- 
sults. In Fig. Q] we show the allowed regions from 
the combined analysis of solar and KamLAND 
data. 
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Figure 2. sin 2 #atm~A to atm regions allowed at 
90%, 95%, 99%, and 3cr C.L. for 2 d.o.f. (the un- 
shaded regions include atmospheric data only). 
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Figure 1. Regions allowed by solar+KamLAND 
data at 90%, 95%, 99%, and 3<r C.L. for 2 d.o.f.. 
Unshaded regions correspond to solar data only. 



2.2. Atmospheric + K2K 

The zenith angle dependence of the /z-like at- 
mospheric neutrino data from the Super-K ex- 
periment provided the first evidence for neutrino 
oscillations in 1998, an effect confirmed also by 
other atmospheric neutrino experiments [2J. The 
dip in the L/E distribution of the atmospheric 
survival probability seen in Super-K gives a 



clearer signature for neutrino oscillations. 

The analysis summarized below includes the 
most recent charged-current atmospheric neu- 
trino data from Super-K, with the e-like and fi- 
like data samples of sub- and multi-GeV con- 
tained events grouped into 10 zenith-angle bins, 
with 5 angular bins of stopping muons and 10 
through-going bins of up-going muon events. No 
information on v T appearance, multi-ring \i and 
neutral-current events is used since an efficient 
Monte-Carlo simulation of these data would re- 
quire more details of the Super-K experiment, in 
particular of the way the neutral-current signal 
is extracted from the data (see Refs. |19l23j for 
details). Here the new three-dimensional atmo- 
spheric neutrino fluxes given in |11| are used, 
in contrast to previous analyses using the Bar- 
tol fluxes [21] ■ With this one obtains the re- 
gions of two-flavour — > v r oscillation parame- 
ters sin 2 #atm and Ato atm shown by the hollow 
contours in Fig. |5J One notes that the Am 2 ATM 
values obtained with the three-dimensional atmo- 
spheric neutrino fluxes are lower than obtained 
previously JH], i n excellent agreement with the 
results of the Super-K collaboration . 

The KEK to Kamioka (K2K) long-baseline 
neutrino oscillation experiment A. probes dis- 
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appearance in the same Am 2 region as probed 
with atmospheric neutrinos. The neutrino beam 
is produced by a 12 GeV proton beam from the 
KEK proton synchrotron, and consists of 98% 
muon neutrinos with a mean energy of 1.3 GeV. 
The beam is controlled by a near detector 300 m 
away from the proton target. Information on neu- 
trino oscillations is obtained by the comparing 
this near detector data with the content of the 
beam observed by the Super-Kamiokande detec- 
tor at a distance of 250 km. 

The data K2K-I sample has been collected in 
the period from June 1999 to July 2001 (4.8 x 10 19 
protons on target) and gave 56 events in Super- 
K, whereas 80.1^5' 2 were expected for no oscil- 
lations. The probability that the observed flux 
is explained by a statistical fluctuation without 
neutrino oscillations is less than 1% 0]. K2K- 
II started in Fall 2002, and released data at 
the Neutrino2004 conference jlj corresponding to 
4.1 x 10 19 protons on target, comparable to the 
K2K-I sample. Altogether K2K-I and K2K-II 
give 108 events in Super-K, to be compared with 
150.91 

lo.o expected for no oscillations. Out of 
the 108 events 56 are so-called single-ring muon 
events. This data sample contains mainly muon 
events from the quasi-elastic scattering -t- p — > 
H + n, and the reconstructed energy is closely 
related to the true neutrino energy. The K2K 
collaboration finds that the observed spectrum is 
consistent with the one expected for no oscilla- 
tion only at a probability of 0.11%, whereas the 
spectrum predicted by the best fit oscillation pa- 
rameters has a probability of 52% 

The re-analysis of K2K data given in ^2] uses 
the energy spectrum of the 56 single-ring muon 
events from K2K-I + K2K-II (unfortunately not 
the full K2K data sample of 108 events, for lack 
of information outside the K2K collaboration). 
Under reasonable assumptions, one can fit the 
data divided into 15 bins in reconstructed neu- 
trino energy. One finds that the neutrino mass- 
squared difference indicated by the i/ M disappear- 
ance in K2K agrees with atmospheric neutrino 
results, providing the first confirmation of oscil- 
lations with Am^ TM from a man-made neutrino 
source. However K2K gives a rather weak con- 
straint on the mixing angle due to low statistics 



in the current data sample. 

The shaded regions in Fig. are the allowed 
(sin 2 (9 atMi ^ m ATM) regions that follow from 
the combined analysis of K2K and Super-K at- 
mospheric neutrino data. One sees that, al- 
though the determination of sin $atm is com- 
pletely dominated by atmospheric data, the K2K 
data start already to constrain the allowed region 
of Am^ TM . Note also that, despite the downward 
shift of Am^rpjyj implied by the new atmospheric 
fluxes, the new result is statistically compatible 
both with the previous one in |19j and with the 
value obtained by the Super-K L/E analysis |2j. 
Note that the K2K constraint on Am\ TM from 
below is important for future long-baseline ex- 
periments, as such experiments are drastically af- 
fected if Am^ TM lies in the lower part of the 3cr 
range indicated by the atmospheric data alone. 

3. THREE-NEUTRINO OSCILLATIONS 

The effective leptonic mixing matrix in gauge 
theories of massive neutrinos was systematically 
studied in 6 . For high-scale seesaw models this 
matrix is approximately unitary and, for three 
neutrinos, may be taken as |26| 3 

K = ^23^13^12 (1) 

where each factor, for example, 

/ ci 3 e^sia \ 
W13 = 1 

\ -e-^sis cis / 

contains an angle and a CP phase. Two of 
the three angles are involved in solar and atmo- 
spheric oscillations, so we set 6 12 = #sol and 
623 = Oatm- All these three phases are phys- 
ical one corresponds to the one present in 
the quark sector (Dirac-phase) and affects neu- 
trino oscillations, while the other two are asso- 
ciated to the Majorana nature of neutrinos and 
show up in neutrinoless double beta decay and 
other lepton-number violating processes, but not 
in conventional neutrino oscillations |32I33| . 

3 This also holds in radiative models |27I28| and models 
where supersymmetry is the origin of neutrino mass 1291 . 
In inverse seesaw models I3UI deviations from unitarity 
may be phenomcnologically important 1311 . 
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parameter 


best fit 


3<T range 


Am^ [lO' 5 


eV 2 ] 


8.1 


7.2-9.1 


Am§ x [1CT 3 


cV 2 ] 


2.2 


1.4-3.3 


sin 2 




0.30 


0.23-0.38 


sin 2 6*23 




0.50 


0.34-0.68 


sin 2 0i3 




0.000 


< 0.047 



Table 1 

Three-neutrino oscillation parameters from the 
global data analysis given in |14j . 



Current neutrino oscillation experiments are 
insensitive to CP violation, thus all phases will 
be neglected (future neutrino factories aim at 
probing the effects of the Dirac phase |34| ) . 
In this approximation three-neutrino oscilla- 
tions depend on the three mixing parameters 
sin 2 0i2, sin 2 023, sin 2 0i3 and on the two mass- 
squared differences Amg OL = Am|i = m\ — to 2 
and Am^rpj^ = Am|i = m| — m\ characterizing 
solar and atmospheric neutrinos. The hierarchy 
ArrigQ L <C Arn^ TM implies that one can set, to a 
good approximation, Amg OL = in the analysis 
of atmospheric and K2K data, and Am^ TM to 
infinity in the analysis of solar and KamLAND 
data. The relevant neutrino oscillation data in 
a global three-neutrino analysis are those of sec- 
tions 12.11 and 12.21 together with the constraints 
from the CHOOZ and Palo Verde reactor exper- 
iments 3 "> . 

The results of the global three-neutrino 
analysis are summarized in Fig. and in 
Tab. n In the upper panels of the figure the 
Ax 2 is shown as a function of the parame- 
ters sin 2 0i2, sin 2 023, sin 2 0i3, Am^i, Am^, min- 
imized with respect to the undisplayed param- 
eters. The lower panels show two-dimensional 
projections of the allowed regions in the five- 
dimensional parameter space. The best fit values 
and the allowed 3a ranges of the oscillation pa- 
rameters from the global data are summarized in 
Tab.Q] This table gives the current status of the 
three-flavour neutrino oscillation parameters. 

In a three-neutrino scheme CP violation disap- 
pears when two neutrinos become degenerate (i 
or when one angle vanihes, 0i3 — ► |36| . Genuine 




{sin 2 e 12 , sin 2 e 33 } s ' n20 i3 



Figure 3. Three-neutrino regions allowed by the 
world's oscillation data at 90%, 95%, 99%, and 3er 
C.L. for 2 d.o.f.Top panels give A% 2 minimized 
with respect to all undisplayed parameters. 



three-flavour effects are associated to the mass hi- 
erarchy parameter a = Arng OL / Am™ and the 
mixing angle 0i3. The left panel in Fig. 0] gives 
the parameter a as determined from the global 
X 2 analysis of 14 . The figure also gives A^ 2 as a 
function of the parameter combination a sin 20i2 
which, to leading order, determines the long base- 
line v e — > Vu oscillation probability |37I38| . The 
last unknown angle in the three-neutrino leptonic 
mixing matrix is 0i3, for which only an upper 
bound exists. The leftt panel in Fig. ^ gives 
Ax 2 as a function of sin 2 0i3 for different data 
sample choices. One finds that the new data 
from KamLAND have a surprisingly strong im- 
pact on this bound. Before the KamLAND-2004 
data the bound on sin 2 0i3 from global data was 
dominated by the CHOOZ reactor experiment, 
together with the determination of Am^ from at- 
mospheric data. However, with the KamLAND- 
2004 data the bound becomes comparable to the 
reactor bound, and contribute significantly to the 
final global bound 0.022 (0.047) at 90% C.L. (3<r) 
for 1 d.o.f. This improved sin 2 0i3 bound follows 
from the strong spectral distortion found in the 
2004 sample 11. Note that, since the CHOOZ 
bound on sin 2 0i3 deteriorates quickly as Am^ TM 
decreases (see Fig. |SJ, the improvement is espe- 
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Figure 4. Determination of a = Arng OL /Am^ TM 
and bound on sin 2 # 13 from current neutrino os- 
cillation data. 



cially important for lower Am^ TM values, as im- 
plied by the new three-dimensional atmospheric 
fluxes . In Fig. [3] we show the upper bound 
on sin 2 #13 as a function of Am^ TM from CHOOZ 
data alone compared to the bound from an analy- 
sis including solar and reactor neutrino data. One 
finds that, although for larger Am^ TM values the 
bound on sin 2 #13 is dominated by CHOOZ, for 
Am ATM £s 2 x 10~ 3 cV 2 the solar + KamLAND 
data start being relevant. The bounds implied by 
the 2002 and 2004 KamLAND data are compared 
in Fig. In addition to reactor and accelerator 
neutrino oscillation searches (Lindners' talk), fu- 
ture studies of the day/night effect in large water 
Cerenkov solar neutrino experiments like UNO 
or Hyper-K |39j may improve the sensitivity on 

sin 2 6>i 3 ma. 

4. ABSOLUTE NEUTRINO MASSES 

Neutrino oscillation data are sensitive only to 
mass differences, not to the absolute neutrino 
masses. Nor do they have any bearing on the 
fundamental issue of whether neutrinos are Dirac 
or Majorana particles |32I33| . The significance 
of the (3[3ol> decay is given by the fact that, in a 
gauge theory, irrespective of the mechanism that 
induces 00ov , it is bound to also yield a Ma- 
jorana neutrino mass 41., as illustrated in Fig. 
El Quantitative implications of the "black-box" 
argument are mo del- dependent, but the theorem 
holds in any "natural" gauge theory. 



Figure 5. Upper bound on sin 2 #13 (1 d.o.f.) 
from solar+KamLAND+CHOOZ data versus 
Am^ TM . Dashed (solid) curves correspond to 
90% (3er) C.L. bounds, thick curves include 
KamLAND-2004 data, thin ones do not. Light 
(dark) regions are excluded by CHOOZ at 90% 
(3cr) C.L. The horizontal line corresponds to the 
current Ato 2 ^^ best fit value, hatched regions 
are excluded by atmospheric + K2K data at 3u. 



Now that oscillations are experimentally con- 
firmed we know that must be induced by 
the exchange of light Majorana neutrinos. The 
corresponding amplitude is sensitive both to the 
absolute scale of neutrino mass as well as the two 
Majorana CP phases that characterize the mini- 
mal 3- neutrino mixing matrix [H] . Fig. shows 
the estimated average mass parameter character- 
izing the neutrino exchange contribution to 
versus the lightest neutrino mass. The upper 
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Figure 6. Equivalence between f3(3ov and Majo- 
rana mass in gauge theories |41| . 
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Figure 7. Neutrino-mass-induced PPov from cur- 
rent oscillation data versus current and projected 
experimental sensitivities. 



(lower) panel corresponds to the cases of normal 
(inverted) neutrino mass spectra. The calculation 
takes into account the current neutrino oscillation 
parameters from |14j and the nuclear matriz ele- 
ments of |16| . In contrast to the normal hierarchy, 
where a destructive interference of neutrino am- 
plitudes is possible, the inverted neutrino mass 
hierarchy implies a "lower" bound for the flPov 
amplitude. Quasi-degenerate neutrinos |42l43j 
such as predicted in |T7|. give the largest (3(3q v 
amplitude. Future experiments |44j will provide 
an independent confirmation of the Heidelberg- 
Moscow results 03j and extend the sensitivity to 
inverse hierarchy models. 

Complementary information on the absolute 
scale of neutrino mass comes from beta decays 
searches as well as cosmology |Tf| . 

5. FOUR-NEUTRINO OSCILLATIONS? 

In addition to the strong evidence for oscilla- 
tions of solar and atmospheric neutrinos there is 
also a hint for oscillations from LSND This 
accelerator experiment at Los Alamos observed 
87.9 ± 22.4 ± 6.0 excess events in the — ► D e 
appearance channel, corresponding to a transi- 
tion probability of P = (0.264 ± 0.067 ± 0.045)%, 
which is ~ 3.3 a away from zero. To explain 



this signal with neutrino oscillations requires a 
mass-squared difference Ato^ snd ~ 1 eV 2 , a 
value inconsistent with the solar+KamLAND and 
atmospheric+K2K experiments described above. 
Four-neutrino schemes 42 49 50 where a sterile 
neutrino is added to the three active ones might 
provide the additional mass scale needed to rec- 
oncile the LSND evidence. 

The updated analysis described in ^31 includes, 
in addition to LSND, solar+KamLAND and at- 
mospheric+K2K data, also the data from short- 
baseline (SBL) accelerator and reactor experi- 
ments with no evidence for oscillations. The os- 
cillation data are divided into the four sets X = 
SOL, ATM, NEV, LSND and the PG method is 
used to test their statistical compatibility in a 
given mass scheme. Following Ref. |51| we con- 
sider the contributions to Axx = Xx ~ (Xx)min 
from different data samples (see |52j for more de- 
tails). In Tab. |21 we show the contributions of the 

4 data sets to Xpg = Xmin for ( 3+1 ) and ( 2 + 2 ) 
oscillation schemes. As expected we note that in 
(3+1) schemes the main contribution comes from 
SBL data due to the tension between LSND and 
NEV data in these schemes. For (2+2) schemes 
a large part of Xpg comes horn solar and at- 
mospheric data, due to the rejection against a 
sterile neutrino contribution of these two data 
sets. The contribution from NEV data in (2+2) 
comes mainly from the tension between LSND 
and KARMEN [53], which does not depend on 
the mass scheme. The parameter goodness of fit 
is obtained by evaluating Xpg f° r 4 d.o.f. |B"4] . 
Using an improved goodness of fit method espe- 
cially sensitive to the combination of data sets 
one finds that (2+2) schemes are ruled out at the 
4.9cr level. The inclusion of MACRO data in the 
analysis would further improve the degree of re- 
jection of these schemes. Note that such a strong 
rejection of (2+2) schemes holds irrespective of 
whether LSND is confirmed or not. 

Although none of the four-neutrino schemes 
provides a good fit to the global oscillation data 
including LSND, it is interesting to consider the 
relative status of the three hypotheses (3+1), 
(2+2) and the three-active neutrino scenario 
(3+0). This is done by comparing the x 2 value 
of the best fit point - which occurs for the (3+1) 





SOL ATM LSND NEV 


Xpg 


PG 


(3+1) 
(2+2) 


0.0 0.4 5.7 10.9 
5.3 20.8 0.6 7.3 


17.0 
33.9 


1.9 x 10" 3 (3.1a) 
7.8 x 10" 7 (4.9<r) 



Table 2 

Contributions of different data sets to Xpg m (3+1) an d (2+2) neutrino mass schemes. 
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Figure 8. Upper bound on sin 2 26>lsnd from 
NEV, atmospheric and K2K data in (3+1) 
schemes. The dotted line corresponds to the 
99% C.L. bound without K2K and using the 1-d 
Bartol flux. Also shown are the regions allowed 
at 99% C.L. (2 d.o.f.) from global LSND data 
and decay-at-rest (DAR) LSND data. 



incompatible with the signal observed in LSND at 
the 95% C.L. Marginal overlap regions between 
the bound and global LSND data exist only if 
both are taken at 99% C.L. Restricting to the 
decay-at-rest LSND data sample jlSj makes the 
disagreement even more severe. This shows that 
4-neutrino descriptions of LSND do not provide a 
satisfactory fit to the world's neutrino data. Al- 
though a reasonable 5-neutrino fit is possible, the 
confirmation of LSND is essential. 



scheme - with the ones corresponding to (2+2) 
and (3+0). First we observe that (2+2) schemes 
are strongly disfavored with respect to (3+1) with 
a A\ 2 = 16.9. For 4 d.o.f. this is equivalent to 
an exclusion at 3.1er. Furthermore, one finds that 
(3+0) is disfavored with a A% 2 = 17.5 (corre- 
sponding to 3.2cr for 4 d.o.f.) with respect to 
(3+1), reflecting the high statistical significance 
of the LSND result. 

On the other hand (3+1) spectra are disfavored 
by the disagreement of LSND with short-baseline 
disappearance data of CDHS and Bugey, leading 
to a marginal GOF of 1.9 x 10~ 3 (3.1cr). Should 
LSND be confirmed a positive signal is predicted 
right at the edge of the current sensitivity. In 
Fig. |S] we give the upper bound on the LSND 
oscillation amplitude sin 2 2#lsnd f rom the com- 
bined analysis of NEV and atmospheric neutrino 
data. From the figure one sees that the bound is 



6. BEYOND OSCILLATIONS 

Non-standard physics may in principle affect 
neutrino propagation properties and detection 
cross sections Such interactions (NSI, for 

short) are a natural outcome of many neutrino 
mass models 56 and can be of two types: flavour- 
changing (FC) and non- universal (NU). They 
may arise from a nontrivial structure of charged 
and neutral current weak interactions with non- 
unitary lepton mixing matrix as expected in 
inverse seesaw models |30I31| . in radiative neu- 
trino mass models |27l28j and in supersymmetric 
unified models |S7j. NSI may be schematically 
represented as effective dimension-6 terms of the 
sub- weak strength sGf, see Fig- El 

In the presence of NSI, the Hamiltonian de- 
scribing atmospheric neutrino propagation has, in 
addition to the standard oscillation part, another 
term -Hnsi accounting for an effective potential 
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Figure 10. Deterioration of 6*13 at a neutrino 
factory in the presence of NSI, explanation in |61| 



induced by the NSI with matter: 



ff NS i = ±V2G F N 



£ 



(2) 



Here +(— ) holds for neutrinos (anti-neutrinos) 
and e and s' parameterize the NSI: \/2GpNfe 
is the forward scattering amplitude for the FC 
process + f — > v T + f and \f2GpNfe' repre- 
sents the difference between + / and jy r + / 
elastic forward scattering. Here AT/ is the num- 
ber density of the fermion / along the neutrino 
path. 

The impact of non-standard neutrino interac- 
tions on the determination of atmospheric neu- 
trino parameters Am^ TM and sin 2 #atm was con- 
sidered in Ref. jSSj treating the NSI strengths 
as free phenomenological parameters and tak- 
ing for / the down-type quark, for definiteness. 
This analysis takes into account both the effect 
of Vy, — > v r oscillations as well as the existence 
of non-standard neutrino-matter interactions in 
this channel. It is shown that, in the 2-neutrino 
approximation, the determination of the atmo- 
spheric neutrino oscillation parameters Am^ TM 
and sin 2 26>atm is practically unaffected by the 
presence of NSI. Future neutrino factories will 
substantially improve this bound |59| . In con- 
trast, as shown in |60| the determination of so- 
lar neutrino parameters is still not yet fully ro- 
bust against the presence of NSI. Before conclud- 
ing, let me mention that, even a small residual 
non-standard interaction of neutrinos in the e — r 
channel leads to a drastic loss in sensitivity in the 
determination #13 through the so-called "golden 
channels" at a neutrino factory |ST], see Fig. I1UI 



This may only be partly overcome by combining 
if different baselines. Only by rejecting NSI at a 
near detector can one improve sensitivities in #13. 

7. CONCLUSION 

Experiment is ahead of theory in neutrino 
physics: despite the great progress achieved re- 
cently we are very far from a "road map" to the 
ultimate theory of neutrino properties. We have 
no idea of the underlying neutrino mass genera- 
tion mechanism, its characteristic scale or its fla- 
vor structure. We have still a long way to go! 

REFERENCES 



7 



T. Mitsui, these proceedings. 
T. Kajita, these proceedings. 
G. Gratta, Neutrino 2004 proceedings. 
T. Nakaya, Neutrino 2004 proceedings. 
S. Weinberg, Phys. Rev. D22 1694 (1980). 
J. Schechter and J. W. F. Valle, Phys. Rev. 
D22, 2227 (1980); D25, 774 (1982). 
M. Fukugita and T. Yanagida, Phys. Lett. B 
174 (1986) 45. 

8. L. Wolfenstein, Phys. Lett. B 107 77 (1981); 
J. Schechter and J. W. F. Valle, Phys. Rev. 
D24, 1883 (1981), Err. D25, 283 (1982); see 
also talks by L. Wolfenstein and J. W. F. Valle 
at Joe's Fest, AIP Conf. Proc. 687, 16 (2003) 

9. M. Goodman, Neutrino Oscillation Industry 
Web-Page, http://neutrinooscillation.org/ 

10. J. N. Bahcall and M. H. Pinsonneault, Phys. 
Rev. Lett. 93, 121301 (2004) 

11. M. Honda, T. Kajita, K. Kasaharaand S. Mi- 



dorikawa, astro-ph/0404457 

12. S. P. Mikheev and A. Y. Smirnov, Sov. J. 
Nucl. Phys. 42, 913 (1985). 

13. L. Wolfenstein, Phys. Rev. D17, 2369 (1978). 

14. M. Maltoni, T. Schwetz, M. A. Tortola and 
J. W. F. Valle, New J. Phys. 6 (2004) 122 
hep-ph/0405172 and references therein. 

15. S. Goswami, talk at Neutrino 2004. 

16. S. M. Bilenky, A. Faessler and F. Simkovic, 
Phys. Rev. D70 033003 (2004) 

17. K. Babu, E. Ma and J. W. F. Valle, Phys. 
Lett. B552 207 (2003); M. Hirsch et al, Phys. 
Rev. D 69 093006 (2004) 



10 



18. 



19. 

20. 
21. 

22. 

23. 

24. 

25. 
26. 

27. 
28. 
29. 

30. 

31. 



C. H. Albright, hep-ph/0407155; W. Grimus 
et al, piep^p h/0408123| S~~M. Barr and 
I. Dorsner, Nucl. Phys. B 585 79 (2000); 
H. S. Goh et al, Phys. Lett. B 587 105 (2004); 
M. Raidal, |hep-ph/0404046| H. M inakata and 
A. Y. Smirnov, |hep^ph /0405088 , G. Altarelli 
and F. Feruglio, Proc. of 10th Int. Workshop 
on Neutrino Telescopes, Venice (2003) 237 
M. Maltoni et al, Phys. Rev. D68, 113010 
(2003); Phys. Rev. D67, 013011 (2003) 
G. Fogli et al, Phys. Rev. D66, 053010 (2002) 
A. B. Balantekin and H. Yuksel, Phys. Rev. 
D68, 113002 (2003) 

G. Fiorentini et al, Phys. Lett. B558, 
15 (2003); H. Nunokawa, W. C. Teves and 
R. Zukanovich, JHEP 11, 020 (2003) 
M. Gonzalez-Garcia et al, Phys. Rev. D63, 
033005 (2001) 

G. Barr, T. K. Gaisser and T. Stanev, Phys. 
Rev. D39, 3532 (1989). 



32. 



33. 



34. 



35. 



3G. 

37. 
38. 



http: / / eps2003.physik.rwth-aachen.de 
Data Group, 



al, 



Particle Data Group, K. Hagiwara et 
Phys. Rev. D66, 010001 (2002). 
A. Zee, Phys. Lett. B93, 389 (1980). 
K. S. Babu, Phys. Lett. B203, 132 (1988). 
M. Hirsch and J. W. F. Valle, New J. Phys. 
6, 76 (2004) |hep-ph/0 405015 
R. N. Mohapatra and J. W. F. Valle, Phys. 
Rev. D34, 1642 (1986). 

J. Bernabeu et al, Phys. Lett. B187, 303 
(1987); G. C. Branco, M. Rebelo and J. W. F. 
Valle, Phys. Lett. B225, 385 (1989); N. Rius 
and J. W. F. Valle, Phys. Lett. B246, 
249 (1990); F. Deppisch and J. W. F. Valle, 
hep-ph/0406040. 

J. Schechter and J. W. F. Valle, Phys. Rev. 
D23, 1666 (1981). 

M. Doi et al, Phys. Lett. B102, 323 (1981); 
S. Bilenky et al, Phys. Lett. B 94 49 (1980) 
K. Dick, M. Freund, M. Lindner and A. Ro- 
manino, Nucl. Phys. B562, 29 (1999) 
M. Apollonio et al, Phys. Lett. B466, 415 
(1999); F. Boehm et al, Phys. Rev. D64, 
112001 (2001) 

J. Schechter and J. W. F. Valle, Phys. Rev. 
D21, 309 (1980). 

M. Freund, Phys. Rev. D64, 053003 (2001) 
E. K. Akhmedov et al, hep-ph/0402"T75| 



39. 



40. 



41. 

42. 

43. 

44. 
45. 
46. 
47. 



C. Yanagisawa, Proc. of AHEP2003, pub- 
lished at JHEP, PRHEP-AHEP2003/062, ac- 



49. 



50. 



51. 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



61. 



Mohapatra, Phys. 
Valle, Phys. Lett. 



cessible from http://ific.uv.es/ahep/ 
E. K. Akhmedov et al, JHEP 0405 057 
(2004); M. Blennow et al, Phys. Rev. D69, 
073006 (2004) 

J. Schechter and J. W. F. Valle, Phys. Rev. 
D25, 2951 (1982) 

D. O. Caldwell and R. N. 
Rev. D48, 3259 (1993). 
A. Ioannisian and J. W. F 

B 332 93 (1994) 

H. V. Klapdor et al, |hep^ph/9910205| 

H. Klapdor etal, Phys. Lett. B 586 (2004) 198 

A. Osipowicz, et al, hep-ex/0109033 

S. Pastor, these volume; S. Hannes- 

tad, |hep-ph/0404239| G. Fogli etal, 

piep-ph/0408045|~ 

A. Aguilar et al, Phys. Rev. D64. 112007 
(2001) 

J. T. Peltoniemi and J. W. F. Valle, Nucl. 

Phys. B406, 409 (1993) 

J. T. Peltoniemi, D. Tommasini and J. W. F. 

Valle, Phys. Lett. B298, 383 (1993). 

M. Maltoni et al, Nucl. Phys. B643, 321 

(2002), |hep-ph/0207T57| . 

M. Maltoni, T. Schwetz and J. W. F. Valle, 

Phys. Rev. D65, 093004 (2002) 

E. D. Church, K. Eitel, G. B. Mills and 
M. Steidl, Phys. Rev. D66, 013001 (2002) 
M. Maltoni and T. Schwetz, Phys. Rev. D68, 
033020 (2003) 

S. Pakvasa and J. W. F. Valle, Proc. of the 

Indian Nat. Acad, o f Sciences, Vol. 70 A, No.l, 

p.189 - 222 (2004), |hep-ph/030106ll 

J. W. F. Valle, Prog. Part. Nucl. Phys. 26, 

91 (1991); for a good example see Ref. 

L. J. Hall, V. A. Kostelecky and S. Raby, 

Nucl. Phys. B267, 415 (1986). 

N. Fornengo et al, Phys. Rev. D65, 013010 

(2002) 

P. Huber and J. W. F. Valle, Ph ys. Lett. 

B523, 151 (2001), hep-ph/0108193 . 

O. G. Miranda, M. A. Tortola and J. W. F. 

Valle, hep-ph/0406280. 

P. Huber, T. Schwetz and J. W. F. Valle, 

Phys. Rev. D66, 013006 (2002); Phys. Rev. 

Lett. 88, 101804 (2002) 



